Transport in Graphene p-n Junctions in Magnetic Field 
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Ballistic transport in graphene p-n junctions in the presence of magnetic field exhibits two distinct 
regimes: At low fields, transport is partially suppressed by the field. When the field exceeds a 
certain critical value, the junction is pinched off by the Landau level formation. Transmission and 
conductance are found in the entire range of fields using Lorentz boost and mapping to the Landau- 
Zener problem. We show that perfect transmission occurs at a field-dependent coUimation angle, 
indicating that the chiral dynamics of massless Dirac fermions persists at a finite magnetic field. A 
current switch, utilizing field-tunable coUimation angle, is proposed. With a generalization of the 
developed approach we study transmission through p-n junctions in graphene bilayer. 



Graphene p-n junctions, fabricated recently in locally 
gated samples 0,11, 01 , provide a new tool to study elec- 
tron transport. Charge carriers in graphene mimic rela- 
tivistic Dirac particles with zero mass and linear disper- 
sion relation e = ±?;i?|p| with vp w lO^cm/s. Graphene 
p-n junctions are predicted to exhibit signatures of chi- 
ral dynamics of massless Dirac particles: perfect trans- 
mission normal to the junction 0, Q and coUimation of 
the transmitted particles . Ballistic transport in p-n 
junctions was proposed as a means to realize an electron 
lens 0. 

The properties of the p-n-p system studied in Ref. [H 
could be mainly explained by conduction in the disor- 
dered p and n regions, rather than in the p-n junctions. 
Likewise, the effects in quantizing magnetic fields [H, Q 
were understood from edge state transport in the p and 
n regions, with the p-n interface merely providingmode 
mixing f§\ . In neither of the experiments [l|, 0,3 the 
effects of ballistic transmission [4, 6„ Q seemed to stand 
out. 

This is not too surprising, given that direct detection 
of the effects ^4, 5, 6, 7J would require an angle-resolved 
measurement and/or very clean samples. Alternatively, 
one can ask if the behavior 0,|3,@i[3| can be inferred from 
the dependence of transport properties on the magnetic 
field that often provides valuable insights into electron 
dynamics. It is interesting therefore to better understand 
the signatures of ballistic transmission in external mag- 
netic field, which is the main purpose of the present work. 

We start by noting that the coupling of an electron 
to external fields reflects relativistic character of charge 
carriers in graphene with the speed of light c replaced by 
vp- In relativistic electro-magnetic theory the fields E 
and B are treated on equal footing, playing the role of 
each other in a moving reference frame. The dynamics 
of a relativistic particle in uniform fields depends only on 
the Lorentz invariants E^— B^, E.B {^j. In particular, the 
dynamics in crossed fields, E.B = 0, can be of two main 
types, magnetic and electric, depending on the relative 
strength of the fields E and B. In the first case, B > 
E, the particle trajectories are described by cyclotron 
motion superimposed with a drift perpendicular to E. In 



the second case, E > B, the trajectories are similar to 
those in the absence of B field, moving asymptotically 
parallel to E and exhibiting no cyclotron motion. 

Quantum transport in these two regimes, magnetic and 
electric, was discussed a while ago 10|, in the context 



of interband tunneling in two-band semiconductor sys- 
tems modeled by the Dirac equation. Naturally, both of 
these regimes can be realized in graphene p-n junctions. 
In the magnetic case, realized for B > {c/vf)E fl^, elec- 
tron motion is described by quantized Landau levels with 
a linear dispersion in the momentum perpendicular to 
E, i.e. parallel to the junction. This defines relativistic 
Quantum Hall edge states transporting charge along 
the p-n interface. Cyclotron frequency in this regime as 
a function of E vanishes at E — {vf/c)B, signaling col- 
lapse of the Landau levels and Quantum Hall effect ^12.]. 

In the electric regime B < {c/vf)E, which will be 
of main interest for us here, electrons can move freely 
along E, transporting electric current through the junc- 
tion. Transmission coefficient, found below as a function 
of S, is shown to vanish at the critical field 



B = B^ = {c/vf)E. 



(1) 



The effect of increasing magnetic field is therefore to 
pinch off transport through the junction, and transform 
it into the edge state transport along the junction in the 
Quantum Hall state at i? > Similar conclusions for 
tunneling suppression by transverse magnetic field in 3D 
junctions modeled by Dir ac p articles with a finite mass 
were obtained in Refs. 

In our approach, we solve the Dirac equation in crossed 
E and B fields exactly with the help of a Lorentz boost. 
This allows us to treat the monolayer and bilayer cases 
on equal footing. We find collimated transmission peaked 
at an angle sitlOb = B/B^,, with unit transmission at the 
peak, 9 = 9b, as in the absence of magnetic field. The 
net conductance, found by integrating transmission over 
angles 9, for a wide junction is given by 



G{B < B,) = 



2TTh d 



(I - {B/B^f) 



2n3/4 



(2) 
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FIG. 1: Schematic of a p-n junction in a locally gated sample. 
For the geometry shown (wide and short sample) the conduc- 
tance is dominated by the junction. Magnetic field suppresses 
conductance as G{B) oc (1 - (B/B.)^)^/*, Eq.©. 



where d = {hvp /\eE\y/^ and w is the p-n interface length 
(see FigH]). The suppression of tunneling ^ precedes 
formation of edge states at the p-n interface at i? > . 

To estimate the critical field B^, for the parameter val- 
ues of Refs.P, [2, 3] one would have to account for screen- 
ing of the in- plane field created by gates [3] . To bypass 
these complications, we assume that a density variation 
of order uq ^ 10^^ cm~^ is created in a p-n junction 
across a distance £ « 50 nm. Then the field felt by the 
electrons is eE ~ hvp^T^no/i, giving 



= {c/vf)E - {hc/e)^ 



(3) 



In terms of the magnetic length is = \J hc/eB, this 
translates into l\ = t/ ^J-ktlq « 260 nm^, yielding an 
experimentally convenient value of ^ 2.5 T. 

For the p-n junction contribution to dominate over the 
conduction in the p and n regions, it is beneficial to be in 
the ballistic regime, similar to Refs.[l, 0, Q, and to use 
wide and short samples (see Fig[T|). These requirements 
are more relaxed for p-n junctions in epitaxial and bilayer 
systems, where tunneling is exponentially suppressed ow- 
ing to the presence of a spectral gap (see below). 

We first consider transport in the p-n junction in the 
absence of magnetic field. Massless Dirac particles in 
graphene moving near the p-n interface in a uniform in- 
plane electric field are described by the Hamiltonian 



H = e<^(x) + vf£, 



p+ 
n_ 



p±^pi±ip2, (4) 



where is the electrostatic potential used to create the 
junction, and ^ = ±1 for the points K and K' . We 
consider a p-n interface parallel to the x axis (Fig[T]), 
with the external field E || y described by <y5(x) = —Ex2- 
The eigenstates of (jj]) are characterized by the mo- 
mentum component parallel to the junction, ip{t,'x.) = 
g-iet+ipia:i^^^^^^ giving a ID problem for i){x2)- Follow- 
ing Ref . (isl ■ we choose to write this problem in momen- 
tum representation 

- ieE dtp/dp2 = H^, H = vpipicri - P2cr2) - £■ (5) 

As noted in Ref. [l5| , momentum representation provides 
direct access to the asymptotic scattering states, and is 
thus more beneficial than the position representation. 



Indeed, Eq.([5]), interpreted as a time-dependent evolu- 
tion with the Hamiltonian H, "time" p2, and "Planck's 
constant" eE, can be identified with the Landau-Zener 
problem for a two-level system evolving through an 
avoided crossing. Hence the probability to be transmit- 
ted (reflected) in the Dirac problem translates into the 
probability of a diabatic (adiabatic) Landau-Zener tran- 
sition. The transmission coefficient can thus be found 
using the answer for the latter [l3| , giving 



T{pi) = exp(-7rft?;FPi/|e£^|), 



(6) 



which agrees with the results of [1, [l^ (see also [l6|). 

Alternatively, the result (O can be put in the context of 
Klein tunneling that links transmission of a Dirac particle 
through a steep barrier with electron/hole pair creation. 
The pair creation rate can be found as the probability of 
an interband transition occuring when the particle mo- 
mentum evolves as p2 = eEt. Because each created pair 
transfers one electron charge across the p-n interface, the 
pair creation rate is equal to the tunneling current. 

To analyze transport in the p-n junction in the pres- 
ence of a magnetic field, it will be convenient to rewrite 
the Dirac equation ([4]) in a Lorentz-invariant form 



(7) 



where 7^ are Dirac gamma-matrices, 7^^ — (73, 
= —ia2, 7^ = — icTi, and ■0 is a two-component 
wave function. Here we use the space-time nota- 
tion for coordinates = {vpt, xi,X2), momenta 
Pfi = f>'{iv~j^^dt,—idxi,—idx2)i and external field — 
(ao, fli, a2). The fields E || y and B || z are described by 



ao 



-Ey, ai 



-By, a2 = 0. 



(8) 



The Dirac equation ([7]) is invariant under the Lorentz 
group {d — 2 + 1): 

xf^'^A^^x^, p^, =a;:,p^, a^,^A';,a^, (9) 
^' = 5(A)0, (10) 

where S{A) = exp (iti;^i,[7^, 7"^]) for ^ = exp(w). 

We first find transmission quasiclasically, using the 
same factorization as above, ip{t,x) = e^^^*'^'^P'-'^'-ip{x2), 
which gives a ID problem for ^p(x2)'. 

{-f°{e + ax)+j\pi + bx)-i-/^d^)ij{x) = (11) 

where a — -^E, b = ^B, x = X2- Eq. pTjl can be cast in 
the form of evolution with a non-hermitian Hamiltonian: 

idxip{x) = ((e + ax)(T2 + i{pi + bx)a3) tpix). (12) 

Now, we apply the adiabatic approximation, constructed 
in terms of cc-dependent eigenstates and eigenvalues of 
the non-hermitian Hamiltonian. The eigenvalues are 
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FIG. 2: Angular dependence of transmission for different 
magnetic field values, T{e) = e-"T'(='"*-^/^*'^ Eq.([T5ll, for 
a = 'jT{d/XF)^ = 20. Transmission reaches unity at a field- 
dependent angle 9b ~ arcsinS/B, . 




FIG. 3: Field-controled switching of coUimated current flow 
through a ballistic p-n junction between different contacts. 



±k{x), where k{x) = (s + ax)^ — {pi + bx)-^. This 
quantity is imaginary in the classically forbidden region 
xi < X < X2, where xi^2 = {e ± pi)/{a ± b). The WKB 
transmission coefficient then equals e^"^, where 



S 



2 / hiiK{x)dx 

'XX 



(piQ - ebf 

(a2 _ 52)3/2 



(13) 



For B = our WKB result ^ agrees with Eq.®. 

The problem (O, ^ can be solved exactly with the 
help of a Lorentz transformation chosen so as to elimi- 
nate the field B. (This is possible because the Lorentz- 
invariant combination B.E equals zero.) For a not too 
large magnetic field, B < B^, = :^E, we can eliminate B 
by a Lorentz boost with velocity parallel to the junction: 



A: 



7 7^ 
7^ 7 
1 



1 



(14) 



Choosing the boost parameter as (3 — —vpB/cE — 
—B/B^,, in the new frame we have B' ^ 0, E' = E/^. 

Because B' = 0, the transmission coefficient for an 
electron with momentum p'l parallel to the p-n junction 
is given by T = g"'^''"^'^ | j^-^ ^^iq new frame (see 
Eq.® and @, 15|). Expressing p'l and E' through the 



quantities in the lab frame, we obtain 



T{pi 



[hvF/\eE\Y'\ (15) 



e = e/vp, which coincides with the WKB result (|T3|) . 

In passing from the moving and lab frames we used 
the fact that the transmission coefficient T, Eq. pS)) . is 
a scalar with respect to Lorentz transformations p4)) . 
This is true because transmission and reflection at the p- 
n interface is interpreted in the same way by all observers 
moving with the velocity parallel to the interface. 

The dependence of the transparency ([TS]) on the elec- 
tric field E is such that T grows as E increases. This 
is a manifestation of the Klein tunneling phenomenon in 
which steeper barriers yield higher transmission. 

The result (fTS]) features exponential suppression of 
tunneling by B field for all momenta except pi = —j3e 
that yields perfect transmission. This corresponds to the 



incidence angle 9b — arcsinS/S* (see Fig 15]). At equal p 
and n densities, the velocities of transmitted particles are 
coUimated at « 0b , with the coUimation angle variance 
determined by Api w (i^^7^'^/^. This gives an estimate 

M^{Xld){l-{B/B,fY'\ X = VF/eF. (16) 

We conclude that the nearly unit transmission, which 
occurs perpendicular to the p-n interface at i? = 4, f|], 
persists at finite magnetic fields, albeit for 6b 0. This 
behavior of the coUimation angle can be used to realize a 
switch (see Fig 13]), in which current is channeled between 
diferent pairs of contacts by varying the B field. 

The p-n junction net conductance can be found from 
the Landauer formula 

where w is the length of the junction interface (see FiglT]), 
and the states contributing to transport are those at the 
Fermi level, e = ep- For a wide junction, w ^ d ^ Xp, 
extending integration over pi to infinity we obtain the 
(1 - [B/B^ff* dependence ©. 

It is interesting to apply these results to epitaxial 
graphene, described by massive Dirac particles e = 
±(w|-p^ + A^V_/^ with an energy gap A induced by the 
substrate ISj, ll9| . The generalization amounts to replac- 



ing byp^ + A^/w|n in Performing Lorentz transfor- 
mation, we find exponential suppression of conductance: 



G{B) 



2 

e w 
2TTh~d 



(1 - P') 



,2n3/4 



exp 



( nd- 

V "1(1- 



^2)1/2 



(18) 



(c/. Refs.[i3,[ii|). The angular dependence of transmis- 
sion in this case is the same as in the massless case. 

We note that G{B > B^,) — does not necessarily 
mean that the system ceases to conduct. The behavior 
predicted by Eq.Q a.t B > B^ should be interpreted as 
2D transport pinching off by the onset of the Quantum 
Hall effect. In that, just the part of the conductance 
proportional to the sample width w vanishes, while the 
edge mode contribution remains nonzero. 
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FIG. 4: Transmission of a bilayer p-n junction vs. momen- 
tum parallel to the interface and vertical field. Note the dou- 
ble hump structure with unit transmission at the peak, and 
suppression of tunneling at large u and pi [E' = O-OSA^/cwf). 



Our approach can be readily generalized to described 
recently fabricated p-n junctions in graphene bilay- 
ers [lOl • The bilayer Hamiltonian 21 1 includes the stan- 
dard monolayer tight-binding part, as well as a direct 
coupling between the adjacent sites B, A of different 
monolayers and a weaker coupling between non-adjacent 
sites A, B: <C 7^^ ~ 0.6 eV in notation of Ref . [2H . 
Here, for simplicity, we ignore 7^^ and denote 7^^ as A. 

It is convenient to write the bilayer Hamiltonian, lin- 
earized near the Dirac points, in pseudospin notation, 
using T3 = ±1 to label the monolayers. The inter- 
layer coupling takes the form TJa = A (t+ct- + T_cr+) = 
^ (rifji -I- T2(T2), where cr± = \{(7i ± 1(12), t± = \{ti ± 
1x2). This gives the Hamiltonian 



H = VFPicri - + ^uTa + y (ticti + T2cr2) , (19) 

where u is the vertical field that opens a gap of size \u\ 
in the bilayer spectrum. Multiplying the time-dependent 
Schrodinger equation by (T3 , we rewrite it as a Dirac equa- 
tion ([7]) with a fictitious r-dependent gauge field: 



g^)ij = 0, 5^ - iir3,-Ari, Ar2 , (20) 



7^(Pp 



where u = u/2vf, A — Is./2vf^ and the external field 
is defined in the same way as above. 

Under Lorentz boost ([T^ the equation ^^{Pfi — — 
g^) = changes covariantly with the momenta and fields 
transforming via p' — Ap, a' — Aa, g' — Ag, giving 
5m' = 2^ {l{uT3 - (3 An), 7(/3Mr3 - An), Arj). Choos- 
ing P so as to eliminate the B field, we find the trans- 
formed Hamiltonian H' — ~eE'x2 + Hk{p'i,p'2), where 



Hk{p[,p'2)^\l{uT3~pAn) 



VFP'l - -jliPuTj, - Ati) ) (71 



VFP2 



1 



(21) 



Ar2 (72 



treat e'lp' = H'lp' as a first-order differential equation 
ieE'dip/dp2 = {Hk{p'i,P2) - e')ip. 

We evaluate the transfer matrix of this equation numer- 
ically, and find that in the physically interesting case 
u <C A, the lowest and the uppermost energy levels of 
Hk are decoupled from the two middle levels. The 4x4 
transfer matrix is thus reduced to a 2 x 2 matrix, yielding 
the transmission and refiection coefficients. 

Transmission features an interesting behavior as a 
function of external fields and particle momentum (see 
FigU]). It has a symmetric double hump profile as a func- 
tion of pi and u vanishing between the humps (unlike 
single gaussian peak in the monolayer case) and, some- 
what unexpectedly, perfect transmission at the peak. At 
large pi and u, because of the energy gap opening, trans- 
mision is strongly suppressed. Conductance, found from 
the Landauer formula p?]) . also exhibits strong suppres- 
sion at increasing u and B, qualitatively similar to the 
gapped monolayer case, Eg. p^ . 

We benefited from useful discussions with C. M. Mar- 
cus and D. A. Abanin. This work is supported by the 
DOE (contract DEAC 02-98 CH 10886), NSF MRSEC 
(DMR 02132802) and NSF-NIRT DMR-0304019. 



Working in the momentum representation, as above, we 
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